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I. INTRODUCTION 

Light propagation and storing in laser driven atomic media have recently become a subject of numerous studies, 
stimulated both by a fundamental character of the discussed problems and by possible applications. A generic system 
for such studies is a medium of three-level atoms interacting with two laser beams (signal and control) in the A 
configuration. The essence of the process of light storing is a coherent mapping of the signal pulse into an atomic 
excitation, described by an atomic coherence, due to a switch-off of a control laser beam. By switching the control 
field again one restores the trapped pulse, preserving all the phase relations. This means in fact first writing down 
an information conveyed by the photons and later reading it back. The whole process can be effectively interpreted 
in terms of a joint medium+field excitation called dark state polariton. For recent reviews on light storage see, e.g., 
Refs 0,EQ. 

Admitting additional fields and thus additional active atomic states enriches the dynamics of the process and gives 
new possibilities of its control. In particular by adding another upper state one obtains a double A system Q in which 
it is possible to make the medium transparent simultaneously for two signal pulses or to split a signal pulse into two 
ones of different frequencies 0, |(| . Channelization of information is also possible in the case of inverted Y systems • 
Coupling the A system with a side level makes it possible to release a pulse of a frequency different from the original 
one or to temporarily prevent the signal from being released Q . Recently the dynamics of light propagation in 5-level 
atomic systems in the M configuration has also been discussed 0. 

A special example is the so-called tripod atomic system in which three low-lying levels are coupled with an upper 
level by a signal field and two control ones. By a manipulation on control fields one can steer the propagation, storing 
and release of the signal, which may be stored in the form of two atomic coherences. Light propagation and slowdown 
in such a medium have been studied in Refs [Tol 1 1 1| while the problem of light storing has been discussed in detail 
in our recent paper [Ti| ■ An attempt of releasing the trapped pulse by an arbitrary combination of two control fields 
results in general in splitting it into two parts one of which was leaving the sample while the other one remained 
trapped and could be released by switching on a new set of control fields. The evolution of such a system could be 
effectively described in terms of a couple of polaritons. 

The previously used classical description of dividing the pulse into two or more parts may be essentially enriched 
by treating the pulse quantum-mechanically and by asking new questions about the photon statistics, especially 
important in the case of nonclassical light states, e.g., Fock states with a small photon number. Such an approach 
is also necessary in the consideration of quantum information storage in an atomic memory. Wang et at |l3l ] 
demonstrated a possibility of time splitting of a single photon by first storing it in one coherence, than by a sequence 
of transferring a part of the excitation into the other coherence and of a photon release from the latter. An additional 
transfer procedure based on a fractional adiabatic population transfer (F-STIRAP) had to be performed before each 
release procedure. 

In this paper we investigate light propagation and storage of a quantum signal field in a medium of atoms in 
the tripod configuration. We discuss mapping of a quantum wave packet into two atomic coherences which means 
in particular a channelization of a single photon. In our approach, being more general than that of Ref. |l3[ . at 
the release stage the two control fields may be twice switched on and off simultaneously. If the control pulses are 
complementary, which means that their amplitude and phase relations are properly chosen, the whole of the stored 
signal is divided into two parts, which are released at two separated time instants. In the case of a single trapped 
photon this leads to its time-entangled state, but without a need of an additional transfer stage of Ref. 13]. We 
discuss in detail the case in which two time-separated photons are trapped at the same place of the sample due to 
an action of two complementary sets of control pulses. They are later released by a different pair of sets of such 
pulses. We examine the possibility of photon coalescence, i.e. a release of both photons at the same release stage. By 
changing localizations and shapes of the stored photon wave pac kets and the details of the control fields we obtain an 
analogue of the Hong-Ou-Mandel interferometer [14L Ha, Ha. Il7| working on stored light. In contradistinction to the 
original version of such an interferometer we have to do with photons propagating always in the same spatial direction 
but well resolved in time. 



II. QUANTUM POLARITONS IN THE TRIPOD CASE 

In order to investigate quantum effects in light storing and propagation both the signal and the medium have to 
be described quantum-mechanically. Some formulae of this section, pertaining to quantum operators, resemble those 
presented in our previous paper |12( for classical pulse. However, they are rederived below not only to make the 
presentation complete but also because they will serve to a discussion of physical effects connected with a quantum 
field, in particular with single photon states. In particular polaritons which in the present formulation arc quantum- 
mechanical operators will play a central role in describing photon correlation effects. 
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Consider a one-dimensional medium of atoms in the tripod configuration (Fig. ^| including an upper state a and 
three stable lower states b (initial), c and d. The sample is irradiated with three collinear laser beams. The quantized 
signal field couples the states b and a and is written as: 

e(z, t) = e (+) (z, t) + e (_) (z, t) = T, k g k a k exp[i(kz - tot)] exp[-i(k ab z - uj ab t)} + h.c, (1) 



where a k is the annihilation operator of a photon with wave number k, g k — <J 2 e„v > ^ * s ^ nc c l uari 'tiza'tion volume, 
eo - the vacuum electric permittivity and k a p — io a p/c = (E a — Ep)/(hc), a,j3 — a,b,c,d. The two control fields 
£2,3 (i) cos(fc2,3Z — u>2,3t) are treated classically and are supposed to be strong enough for the propagation effects to be 
neglected; they couple respectively the states c and d with a. The medium excitation is described in the formalism of 
the second quantization by the flip operators a a p = \a >< /3| exp(ik a pz), k a p being the corresponding wave vector. 
The medium is treated in a continuous way; this is why j- J dz in the above equation has replaced the sum over N 
atoms (L is the length of the sample). 

The interaction hamiltonian in the rotating-wave approximation reads 



N f 

H = — dz(-d ab a ab exp[-i(fc afc z - u ab t)]Y,a k g k exp[i(kz - u)t)] + h. 



c. 



--jdac<Jac exp[-z(fc ac z - w ac t)]e 2 (i) exp[i(fc 2 z - uj 2 t)] + h.c (2) 

-^d ad <Jad exp[-i(k ad z - cu ad t)]e 3 (t) exp[i(k 3 z - w 3 i)] + h.c), 

where d a p are the dipole moment's matrix elements, assumed for simplicity to be real and positive. 

In our analysis in terms of polaritons it is enough to restrict oneself to the deterministic part of the Heisenberg 
evolution (relaxation- and noise- free) . The equations of motion for the flip and signal field operators in the first-order 
approximation with respect to the signal field, in resonance conditions read 

i& ba {z, t) = -^ e (+) (z, t) - Q 2 (t)a bc (z, t) - {l 3 (t)<r bd (z, t), 

i& bc (z, t) = -D,l(t)a ba (z, t), (3) 
ia bd (z,t) = -D,l(t)a ba (z,t), 

(m + c m> {+) ^ t ) = i f9 2 d ba a ba {z,t), 

where g is the value of g k for the central field frequency of the signal field, and SI2 (i) = d ac e2(t) / (2h) and £l 3 (t) = 
d a d^{t) I (2h) are the Rabi frequencies for the control fields. 

It is convenient to express the solutions of Eqs (3) in terms of two field+medium excitations called polaritons: 
^(zjt) (the dark-state polariton) and Z(z,t) 

*(z,i) = exp(-ix){e (+) (z,i) cos6» - ^ sin0[exp(ix 2 ) cos0er bc (z, t) + exp(ix 3 ) sin (j)(J bd {z, t)]}, 

Z(z,t) = j^[exp(ix2) sin (f>a bc (z 7 1) - exp(ix3) cos cj)(T bd (z 7 t)}, (4) 

where the mixing angles 9 and <fi are defined as tan 6* = k/Q, tanc/> = IO3/O2I with \i — ar g(^2), X3 = arg(^3), 
ft = -\/|£l 2 | 2 + |^3| 2 j k 2 = Ng 2 \d ab \^ 2 /h 2 ; x satisfies the equation x = sin 2 9\2- For simplicity we assume that X2 = X3- 
The field and medium components of the polariton VE" are 

e (+) = exp(ix)\I>cos0, (5) 
exp(ix2) cos (j>a bc + exp(ix 3 ) sin (j>a bd = -jgt exp(ix)* sin 9. (6) 

The evolution equations for the two polaritons read 

(^ + c cos 2 6^)* = exp(-ix) sin 9<j)Z, 

jLz = ix 2 Z-exp{ix)sme4>y. (7) 

The operators ^(z,t)/(gy/L) and Z(z,t)/(gy/L) fulfil typical bosonic commutation relations in the first-order ap- 
proximation, in which a bb = 1, a cc = a dd = a al = 0, ([*(*), &{z')\ = [Z{z), Z\z')\ = g 2 L5(z-z r ), [*(«), Z\z')\ = 0), 
so those operators may be considered annihilation operators of the joint field+medium excitation. The essence of the 
approach used in this paper will be treating a medium excitation due to the stored photon as a superposition of two 
excitations, properly suited to the pulse releasing control fields (see Eq. (9) below). 



4 



If the control fields are proportional, i.e. <j> = and xi = \3 = the evolution equations for the polaritons are 
decoupled. The polariton \1/ propagates through the medium without changing its "shape", i.e. the solution is 

#(z,t)=*(z-c / cos 2 6(t')dt',t = 0), (8) 
Jo 

while the polariton Z = const keeps both its shape and localization. 

The above equations constitute a generalization of the approach of Fleischhauer et al. [? ] for a four-level system 
and at the same time of our previous results 12], here formulated in the language of a quantum signal field. Due to 
such an approach we will be able to discuss the effect of the process in terms of photons rather than in the language 
of splitting a classical pulse. 

If at the beginning the control fields corresponding to a given constant mixing angle <jP and given constant \2 ano - 
%3 are strong enough (6 ss 0), the medium is transparent and an incoming photon enters the medium. Switching 
the control field off (9 = ^) results in the medium becoming opaque and the photon becomes trapped in a coherent 
superposition of the atomic coherences. The propagation and storing are described in terms of the polaritons '5° and 
Z°. The former polariton has evolved from being almost purely electromagnetic to become purely atomic while the 
latter one has not taken part in the evolution. The probability amplitudes that the photon has actually been trapped 
in the coherence o~b c and <Jbd are respectively exp(i%§) cos <fi° and exp(ixg) sin <jp . If now the control fields of the same 
mixing angle and the same phases are switched on again the photon will be released with certainty in the same state 
as the initial one. This is mutatis mutandis a repetition of the case of a usual A system discussed before. 

If the control fields to be used at the release stage are characterized by the angles (ft 1 , x\ an d X3 it is proper to 
describe the further evolution in terms of the corresponding new polaritons and Z 1 . The latter pair of polaritons 
is expressed in terms of the former ones as 

4- 1 = {cos ^ cos 4>° exp[z( X 2 - xl)] + sin ^ sin cfP exp[i(xl ~ + 

{cos sin 0° exp[i( X l - X °)] - sin cos <jP exp[i(xl ~ X°)]}^°, (9) 

Z 1 = {sin^ 1 cos <f>°exp[i( X l - *§)] - cos^ 1 sin^ exp[^ - X3)]}* + 
{sin sin 0° exp[i( X l - X°)] + cos 1 cos cjP exp[i( X l - X° 3 )}}Z a . 

The polariton 'J 1 will be released at turn into a photon, while the Z 1 will remain trapped. It follows from the 
form of the polaritons (Eq. (4)) that applying another pair of the control fields with <fr 2 — £ — if) 1 , X \ = x\ + an d 
X% = X3 leads to an exchange of the role of the polaritons: Z 1 becomes VP 2 which will be released and turn into a 
photon. The sets of control fields marked with (1) and(2) can be considered complementary in the sense that their 
subsequent application implies a certain photon release. 

The above results are much more general than those of Ref. |l3j | in which the authors proposed a time splitting 
of a photon by first storing the pulse in one coherence o~b c (here this corresponds to <fP =0), than by pumping a 
part of the excitation into the other coherence o~bd and finally by a pulse release from the latter coherence (here 
(f) 1 = |). The pumping procedure (F-STIRAP) was to be separated from the release stages which required using light 
pulses transverse to both signal and control fields; the latter should be realized in ultracold gases. Here, due to using 
simultaneously pairs of control fields of the same shape pumping and releasing stages coincide. 

III. BEAM SPLITTER AND HONG-OU-MANDEL - TYPE INTERFEROMETER 

The above-described behavior of the atom+field allows one to use the medium and the system of control fields as 
an effective and flexible beam splitter working on stored light. The technical difference between the usual device and 
that of ours is that incoming photons may arrive from only one direction but at one or two different time instants. 
Also the released photons have the same direction but are separated in time. By changing the mixing angle <\> x and the 
phases xi' 1 an< i X3' 1 °f the control fields one can smoothly regulate the amplitudes and phases of the pulse components 
released at the two stages, which would correspond to changing the transmission and reflection rate of a usual beam 
splitter. 

The operation of photon storing can be performed at two stages (input ports). First we trap a first portion of 
incoming photons by switching on and off a pair of control fields characterized by the parameters <t>°,X2 an d xt- After 
the this part of the trapping operation has been finished we may trap a second portion of photons by applying the 
control fields of the mixing angle \ — (ffi, x% + 11 an( i xV)- Note that the second trapping operation does not affect 
the coherences due to the first one. Thus a photon from the first portion is described by the polariton field and 
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that from the second portion - by Z°. A special case is assumed in which <fp = which means that photons have 
been first trapped in the coherence and next - in er^- 

The release operation consists also of two stages (output ports) connected with two pairs of control fields separated 
in time, the first of which is characterized by the mixing angle tfi 1 and \\ an d X3 an( i the second by ^ — <f> , x\ + n 
and xk- The output includes thus two portions of photons, separated in time, and the photon release is complete. 

In the case of trapping and releasing of two photons we have a new realization of the Hong-Ou-Mandcl effect. In 
its original formulation it concerns a 50-50 photon beam splitter with exactly one incoming photon in each of the two 
input ports (incoming photon directions). Due to interference phenomena the probability amplitudes for obtaining 
exactly one photon in each of the two output ports (outgoing photon directions) cancel out. 

In our realization of the input ports single photons, separated in time, are trapped, one in each of the two superpo- 
sitions of the atomic coherences, connected with complementary sets of control fields. In the output ports photons are 
released, again at two stages, from two other superpositions of the two coherences. The analogue of photon incoming 
from two orthogonal directions is their storing at two storing stages. The analogues of reflection and transmission on 
a traditional beam splitter are photon release at the first or second release stage while the parameters characterizing 
the control fields and the atomic system determine the analogues of the reflection and transmission coefficient and 
phase relations. Using the above-mentioned transformations of polaritons at the storing stage and having in mind 
that the incoming polaritons 'J' , Z° (outgoing polaritons "I/ 1 , Z 1 ) become almost purely electric field operators e 1 , e 2 , 
(e 3 , e 4 ) for t — > — oo (t — > oo) we may give the net result for the transformation of the field operators 



-R31 -R32 
R41 R42 / I e / \ £ 4 
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where 



i? 3 i = cos^ 1 cos(/>°exp[i(x2 - X2)] + sin</> 1 sin0° exp[i(xl - x")], 
R 32 = cos^ 1 sin^°exp[i(x2 - X°)] - sint^ 1 cos <j)° exp [i (xl ~ X3)], 
R 41 = sine/) 1 cos</)°cxp[i(x2 - X°)] - cos^ 1 sin0°exp[ix3 ~ X3OL 

R 42 = sin^ 1 sin0°exp[i(x2 - X2)] + co^ 1 cos0°exp[ix3 - x")]- ( n ) 

Let the respective wave packets of the stored excitation be fi(z) and f2{z), their shape being identical with the 
shape of the initial photon wavepackets and their localization depending on the time instants of the switch-off of 
the control fields. The maximum value (unity) of the overlap of the packets s = J f*(z)fk(z)dz corresponds to the 
situation in which two photons with wavepackets of the same shape have been stored exactly at the same place inside 
the sample. 

The field operators corresponding to the input ports arc <3>°(j) = / f*(z)^°(z)dz, Z°(j) = J f* (z)Z°(z)dz. The 
corresponding operators in the output ports are ^(j) = / f* (z)'^ 1 (z)dz, Z x {j) = j f*(z)Z 1 (z)dz. 

The key commutation relations are * ot (j)} = [Z°(j), Z at (j)} = g 2 L, [*°(1), * ot (2)] = [Z°(l), Z ot (2)] = 

g 2 Ls, with analogous relations for VE" 1 and Z 1 . 

The quantum state of the medium after trapping has been accomplished is constructed as due to a creation of two 
excitations due to complementary control fields and characterized by two possibly different wave packets 

IC >= ^* ot (l)Z° + (2)|0 >, (12) 

where the "vacuum" state |0 > means all atoms in the state \b >). The state corresponding to two photons being 
released at the first stage (photon coalescence) is constructed as due to a creation of two excitations corresponding 
to the first release stage characterized by possibly different wave packets 

IC '^?I7TT1JF*" (1) * 1,<2)|0> - (13) 

where a normalizing factor for the two-particle states has been included. Note that the state corresponding to the 
photon coalescence at the second stage is constructed in an analogous way, using the Z operators instead of 

Calculating the projections with the use of the above mentioned commutation relations yields the probability 
amplitude of photon coalescence at the first release stage 

< CilC >= (1 + N 2 ){cos0° COS0 1 expH( X2 ! - xl)] + ^<f>° cx P H(x§ - X3)]} 

{sin^cosc^expHtxa - X2)] ~ cos^ sin^ 1 exp[-i(x§ - xl)]} ( 14 ) 
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If in particular 0° = 0, which means that exactly one photon has been trapped in each coherence, then the photon 
coalescence probability at the first release stage is 

Pcoai(l) = ^(l + |s| 2 )sm 2 2^. (15) 

If the wave packets overlap (s = 1) we obtain for 4* 1 = j (equal amplitudes of the releasing control fields) that 
Pcoai{^) = I- The same value is obtained for the coalescence probability at the second stage. This means that 
the situation is impossible in which exactly one photon is released at each stage. This is analogue of a symmetric 
beam splitter with exactly one photon in each input port. A reduction of the overlap integral leads in a continuous 
way to the situation in which the probability of the latter situation grows to |; this is an analogue of changing the 
length of the arms of the standard HOM interferometer. For example if the wave packets /i^ are Gaussians with 
standard deviations 5\p separated by a distance a, the probability of releasing a single photon at each stage (photon 
noncoalescence) , being our analogue of coincident registration of single photons in the usual realization, is 

p — 1 = \ [1 -4TTir exp ^ 2 (/+^) )] - (16) 

This probability as a function of the distance a exhibits a minimum called the Mandel dip, reaching zero for equal 
packet widths (see Fig. which corresponds to an ideal photon coalescence. On the other hand in the case of 
coinciding packet centers (a — 0) the probability equals (Si — <5 2 ) 2 /[2((5i + <5|)] (see Fig. |3J. Note that it is only for 
equal packets' widths that the coalescence probability drops to zero. 

The effect can also be controlled by choosing the phases of the control field. From Eq. (14) it follows for example 
that for — (f) 1 = j the coalescence probability at the first stage is 

Pcoal(l) = \{l + \s\ 2 ) sin 2 ( X ° - x\ xl + xl). (17) 

Thus the probability of the photons' coalescence strongly depends not only on the shape of the wave packets and the 
relative amplitudes of the control fields but also, in an oscillatory way, on the relative phases of the latter. 

IV. CONCLUSIONS 

We have analyzed the propagation and storing of photons in an atomic medium in the tripod configuration. An 
adiabatic evolution of the system can be described in terms of a couple of polaritons. Using two complementary 
sets of control fields at the storing stage one can stop two photons arriving at different time instants. Applying two 
other complementary sets of the control fields at the release stage one retrieves two time-entangled photons. The field 
operators of the outgoing photons are expressed by those of the incoming ones, depending on the parameters of the 
control fields. We have thus obtained a kind of a beam splitter, operating on a stored light, in which the input and 
output ports correspond to photons arriving and leaving at different time instants. Our analogues of the transmission 
and reflection coefficients as well phase relations can be regulated on demand. We have demonstrated the effect of the 
Hong-Ou-Mandel - type interference and analyzed its dependence on the shapes and positions of the trapped photon 
wavepackets and on the parameters of the control fields. 

Acknowledgments 

The authors thank Mr. Wojciech Wasilewski for stimulating discussions. The work has been supported by Polish 
budget funds allocated for the years 2005-2007 as a research grant No. 1 P03B 010 28. The subject belongs to the 
scientific program of the National Laboratory of AMO Physics in Toruh, Poland. 



7 



M. D. Lukin, Rev. Mod. Phys. 75, 457 (2003). 

A. Andre, M. D. Eisaman, R. L. Walsworth, A. S. Zibrov and M. D. Lukin, J. Phys. B: At. Mol. Opt. Phys. 38, S589 
(2005). 

M. Fleischhauer, A. Imamoglu J. P. Marangos, Rev. Mod. Phys. 77, 633 (2005)). 

M. D. Lukin, PR. Hemmer, and M. O Scully, Adv. At. Mol. Opt. Phys. 42, 347,(2000). 

A.Raczyriski, J. Zaremba, and S. Zielihska-Kaniasty, Phys. Rev. A 69, 043801 (2004). 

Xiong-Jun Liu, Hui Jing, and Mo-Lin Ge, arXiv-quant-ph/040317 (2004). 

A. Joshi and Min Xiao, Phys. Rev. A 71, 041801(R) (2005). 

A. Raczyriski, J. Zaremba, and S. Zieliriska-Kaniasty, Opt Commun. 217, 275 (2003). 
Pengbo Li, Ying Gu, Kai Wang, and Qihuang Gong, Phys. Rev. A 73 032343 (2006). 
E. Paspalakis and P. L. Knight, J. Opt. B: Quantum Semiclass. Opt., 4, S372 (2002). 
I. E. Mazets, arXiv:quant-ph/0407173 vl (2004). 

A. Raczyhski, M. Rzepecka, J. Zaremba, and S. Zielihska-Kaniasty. Opt. Commun. 260, 73 (2006). 
Tun Wang, M. Kostrun, and S. F. Yelin, Phys. Rev. A 70 053822 (2004). 

C. C. Gerry and P. L. Knight, Introductory Quantum Optics, (Cambridge University Press, Cambridge, 2005). 
P. Grangier, Nature 419, 577 (2002). 

C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 95 2044 (1987). 

C. Olindo, M. A. Sagioro, C. H. Monken and S. Padua, Phys. Rev. A 73, 043806 (2006). 

M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 5094 (2000). 



8 



Figures 



FIG. 1: Level and coupling scheme of the tripod system. 



FIG. 2: The probability of photon noncoalescence as a function of the packet's separation a for 5i = 82 solid line, and 5\ = 3^2, 
dashed line. 



FIG. 3: The probability of photon noncoalescence as a function of the ratio 8281 of the packets' widths for a — 0, solid line, 
and a = 2Si , dashed line. 
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